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We investigate the effect of incoherent perturbations on atomic photoionization due to a femtosec- 
ond mid-infrared laser pulse by solving the time-dependent stochastic Schrodinger equation. For a 
weak laser pulse which causes almost no ionization, an addition of a Gaussian white noise to the 
pulse leads to a significantly enhanced ionization probability. Tuning the noise level, a stochastic 
resonance-like curve is observed showing the existence of an optimum noise for a given laser pulse. 
Besides studying the sensitivity of the obtained enhancement curve on the pulse parameters, such 
as the pulse duration and peak amplitude, we suggest that experimentally realizable broadband 
chaotic light can also be used instead of the white noise to observe similar features. The underlying 
enhancement mechanism is analyzed in the frequency-domain by computing a frequency-resolved 
atomic gain profile, as well as in the time-domain by controlling the relative delay between the 
action of the laser pulse and noise. 
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I. INTRODUCTION 

The role of noise on driven quantum systems is a 
subject of considerable interest and importance. Sev- 
eral studies exist in the literature on this broad topic. 
For instance, the noise-induced effects in the nanoscale 
quantum devices such as Josephson junctions [l], 0, the 
macroscopic phase transitions due to quantum fluctua- 
tions [3[, driven multilevel quantum systems under inco- 
herent environment [4( , and stochastic ionization of Ry- 
dberg atoms by microwave noise |5|, to mention just a 
few. In these examples, the action of noise on a system 
can be broadly classified as being of two types: either 
of destructive nature, i.e., noise must be avoided; or of 
constructive nature necessitating its presence. 

It is this nontrivial latter aspect of noise- induced effects 
that has been subject of intense investigation [H, H, H, 0] • 
In particular, the stochastic resonance phenomenon (SR) 
provides a paradigm for the constructive role of noise in 
nonlinear classical as well as quantum systems 0- The 
essence of quantum SR is the existence of an optimum 
amount of noise in a nonlinear system that enhances its 
response to a weak coherent input forcing [H, [^, [kJ EH • 
Despite the diversity of nonlinear dynamics exploiting 
classical SR, most of the quantum mechanical studies of 
the effect have focused on the so-called spin-Boson model, 
which provides an analog of a classical double- well poten- 
tial [Ell, EH EH ■ However, many other physical systems 
exist, particularly atoms or molecules exposed to strong 
laser pulses, where the quantum dynamics can be non- 
linear and therefore added noise could play an important 
role. 

The presence of noise is also worth studying from the 
point of view of steering quantum dynamics of atomic 
or molecular systems [lj, EH, Elj]. Traditionally this is 
achieved with strong, tailored laser pulses by exploiting 
their nonperturbative and nonlinear interaction with the 
atomic systems. In this context, it has been shown that 
weak noise of various origin in multilevel ladder systems 



plays a crucial role [14| . Many scenarios have been dis- 
cussed, such as the need to either cooperate or fight with 
de-coherence in the closed- loop control [lj]], and to en- 
gineer the environment to achieve the steering of quan- 
tum systems towards a desired state [l8| . In the same 
spirit, white shot noise has been used to dissociate di- 
atomic molecules [ijj] . This has implications for the field 
of quantum control. Motivated by the concept of exploit- 
ing noise in nonlinear systems, one can ask the question 
if noise can serve as an extra tool for quantum control. 
Indeed, the concept of the quantum SR effect has not 
been exploited for additional insight in controlling the 
quantum phenomenon. 

In this article, we provide a detailed study of the influ- 
ence of noise in a generic quantum situation, namely the 
photoionization of a single-electron atom interacting with 
an ultrashort laser pulse. We will demonstrate the con- 
ditions under which a resonance-like behavior emerges 
in the stochastic photoionization process. This noise- 
induced phenomenon is studied for a variety of laser 
pulses, from a few optical cycles duration to very long 
ones, and of varying intensities. Furthermore, we sug- 
gest the experimental observability of the effect by em- 
ploying a broadband chaotic light instead of white Gaus- 
sian noise. Lastly, we characterize the underlying gain- 
mechanism in the frequency-domain in order to identify 
the crucial frequency-bands in the broad noise spectrum. 

The article is organized as follows. Section II intro- 
duces our model of the simplest atom interacting with a 
femtosecond laser pulse and white noise, and describes 
our method to solve its stochastic Schrodinger equa- 
tion. In section III, we show the results of the ioniza- 
tion probability (IP) for various combinations of the laser 
pulse and noise. The existence of a stochastic resonance- 
like behavior is quantified using an enhancement fac- 
tor which is computed from IP. The sensitivity of this 
noise-induced effect is tested with laser pulses of vary- 
ing duration and strength, and other types of noise such 
as a realizable broadband chaotic light. To characterize 
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the enhancement mechanism, we compute the frequency- 
resolved gain profile of the driven atom, and study the 
role of relative time-delay between noise and the laser 
pulse. Finally, section IV provides a summary of results 
with our conclusions. 



II. DESCRIPTION OF THE MODEL 



Hydrogen atom interacting with a laser pulse 
and noise 



We consider as an example the simplest single-electron 
atom, i.e., the hydrogen atom. Due to the application of 
an intense linearly polarized laser field F(t), the electron 
dynamics is effectively confined in one-dimension along 
the laser polarization axis [2(| • The Hamiltonian for such 
a simplified description of the hydrogen atonij which is 
here also perturbed by a stochastic force £(t) [2l|, reads 
as (atomic units, H = m = e = 1, are used unless stated 
otherwise) , 

H(x,t) = ^ + V(x)+x{F(t)+Z(t)}, (1) 

where x is the position of the electron and p = —i d/dx 
is the momentum operator. The external perturbations, 
F(t) and £(t), are dipole-coupled to the atom. The po- 
tential is approximated by a non-singular Coulomb-like 
form, 



V(x) = 



1 



(2) 



\J x 2 + a 2 

Such a soft-core potential with parameter a has been rou- 
tinely employed to study atomic dynamics in strong laser 
fields [221 ] . It successfully describes many experimental 
features of multiphoton or tunnel ionization [2(| , and the 
observation of the plateau in higher harmonic generation 
spectra (22[. 

The laser field is a nonresonant mid-infrared (MIR) 
femtosecond pulse described as, 



F(t) = f(t)F sm(ojt + S). 



(3) 



Here Fq defines the peak amplitude of the pulse, u> de- 
notes the angular frequency, and S is the carrier-envelop 
phase. We choose a smooth pulse envelop f(t) of the 
form, 



/(*) 



sin 2 (^/(2r)), 
1, 



* < T 

T<t<T.„ 



cos 2 (7r(i + r - T p )/(2r)), T p -r<t<T p , 

where T p is the pulse duration and r the time for turning 
the field on and off. 

The noise term £(i) is a zero-mean Gaussian white 
noise having the following properties, 
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and noise intensity D [23| . 



2D S(t-t'), 



(4) 
(5) 




0.01 0.02 0.03 0.04 0.05 0.06 
Peak pulse amplitude Fq (a.u.) 

FIG. 1: The ionization probability Pi as a function of the 
laser peak amplitude Fo. Insets show the ionization flux ver- 
sus time for two different pulses (top parts of curves) of am- 
plitudes, Fo = 0.05 and Fo = 0.02, marked by arrows on IP 
curve. Here lo = 0.057, 5 = 0.0, T p = 2Qtx/uj and r = 2tt/w. 



B. Stochastic quantum dynamics 

The presence of the stochastic forcing term in the 
Hamiltonian as described above, makes the quantum evo- 
lution nondeterministic. Thus an averaging over a large 
number of realizations of the stochastic force is required 
in order to produce a statistically meaningful solution 
of the following time-dependent stochastic Schrodinger 
equation, 



Ft 



H(x,t) V{x,t). 



(6) 



For a given realization, the numerical solution of the 
Schrodinger equation amounts to propagating the ini- 
tial wave function Y&q) using the infinitesimal short-time 
stochastic propagator, 



Us(At) = exp -i 



t+At 



H(x,t)dt 



(7) 



One can compute Ut (At) using the split-operator fast 
Fourier algorithm [24| • Details of the method employed 
are described in the Appendix. Successive applications of 
the stochastic propagator U^(At) advance \^q) forward 
in time. 

Note that the initial state |\&o) is always chosen to 
be the ground state of the system having an energy of 
lb = —0.5 a.u.. This is obtained by the imaginary-time 
relaxation method for a? = 2 [2(|. To avoid parasitic 
reflections of the wavefunction from the grid boundary, 
we employ an absorbing boundary [2~jj . 

The ionization flux leaking in the continuum on one 
side, is defined as 12511 . 



Mx R ,t) = Re[y*pV] 3 



(8) 
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where xr is a distant point (typically 500 a.u.) near the 
absorbing boundary. The ionization rate is integrated 
over a sufficiently long time interval to obtain the corre- 
sponding total ionization probability, 



P = 



Jn(x R ,t)dt. 



(9) 



In the following section, we shall use both the ioniza- 
tion flux, and the photoionization yield, to study the in- 
terplay between the laser pulse and noise. From the point 
of view of stochastic resonance-like phenomena, we aim 
at establishing the constructive role of noise in atomic 
photoionization due to a femtosecond laser pulse. 



III. RESULTS 

A. Optimal stochastic enhancement of 
photoionization 

1. Photoionization as a nonlinear effect 

Let us first consider the response of the atom inter- 
acting with a short but strong laser pulse only. Fig. 1 
shows the ionization probability Pi versus the peak pulse 
amplitude Fq for a 20 cycle long MIR laser pulse (ui = 
0.057). This figure shows that with increasing values of 
Fq the ionization probability first increases nonlinearly, 
and then saturates to the maximum value of unity, for 
Fq > 0.05. The behavior of Pi(Fq) is a characteristic 
signature for many atomic and molecular systems inter- 
acting with nonresonant intense laser pulses [22j. 

The laser pulse produces (nonlinear) ionization of 
the atom which is most easily understood, especially 
in the time domain, with the picture of a periodi- 
cally changing tunneling barrier. Ionization flux is pro- 
duced close to those times when the effective potential 
U(x,t) = V(x) + xF(t), is maximally bent down by the 
dipole-coupled laser field. This is illustrated in the inset 
of Fig. 1 with the temporal evolution of the ionization 
flux for laser pulses (shown in the top parts of inset) with 
two different peak amplitudes F = 0.05 and F = 0.02 
(see arrows in Fig. 1). Time- resolved ionization peaks 
separated by the optical period (2n/ui) are clearly vis- 
ible for both peak field amplitudes. In addition, ■/#(£) 
shows a complex interference pattern [inset of Fig. 1] due 
to the modulated Coulomb barrier for Fq = 0.05. How- 
ever, quite strikingly, if Fq is reduced to 0.02 a.u., the 
ionization flux collapses by around five orders of mag- 
nitude as shown in Fig. 1. One can therefore conclude 
that the photoionization dynamics is highly nonlinear, 
and in particular it exhibits a form of "threshold" dy- 
namics where the threshold is created by the condition 
for over-the-barrier ionization. 
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FIG. 2: Ionization flux for a weak laser pulse Fq = 0.02, 
with three values of noise amplitude, (a) yjD = 0.00024, (b) 
0.0015, and (c) 0.018. Background featureless curves (red) 
show the corresponding purely noise-driven (Fo — 0) flux. 
The flux is averaged over 50 realizations. 



2. Ionization induced by noise alone 

Here we look into the possibility of efficiently ionizing 
the atom, when it is subjected to white Gaussian noise 
only. The interaction time of the atom with the noise 
is kept identical to the laser pulse duration T p (see inset 
of Fig. 3). Fig. 2 shows the evolution of the ionization 
flux (Jn(t)) which is averaged over 50 different realiza- 
tions of the noise. One can see that for small noise am- 
plitudes the ionization flux exhibits a featureless curve, 
producing the ionization flux around 10~ 9 . As the noise 
level is increased, the featureless ionization curve rises 
monotonically as shown in Fig. 2(a)-(c). By integrating 
the stochastic ionization flux, one can compute the corre- 
sponding ionization probability P n , which is simply equal 
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FIG. 3: The ionization probability versus the noise amplitude 
\[D for the noise alone P„ (diamonds) and for the laser pulse 
(Fq = 0.02) with noise P;+ n (squares). The points A-C cor- 
responds to the three cases considered in Fig. 2. The dashed 
line shows the limit of vanishing noise, i.e., the probability 
due to the laser pulse alone Pi. Inset: an example of laser 
pulse and noise 



to the area under the curve (Jn(t)). The resulting noise- 
induced ionization probability P n versus the noise ampli- 
tude is shown in Fig. 3. As can be seen the stochastic 
ionization probability rises monotonically with the noise 
level. For ultra-intense noise, such that its strength be- 
comes comparable to the atomic binding field, obviously 
full ionization can be achieved. We should mention that 
similar effects have been observed in other systems, for 
example, the purely noise-induced molecular dissociation 
[l9l ]. and the ionization induced by weak noise of the 
highly excited Rydberg atoms Q. However, in our case 
we consider the atom to be initially in its ground state. 

Although the application of noise alone, or the laser 
pulse alone, can lead to the atomic ionization, we aim 
to study whether a combination of both the laser pulse 
and noise makes the ionization process more efficient as 
compared to the individual cases. 



3. Simultaneous application of the laser pulse and noise 
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FIG. 4: The enhancement in photoionization due to quan- 
tum SR. The points marked A-C correspond to the noise am- 
plitudes of Fig. 2(a)-(c), respectively. At point B the ratio 
\J D opt I Fq = 0.075. The error bars indicate the standard de- 
viation of n calculated using more than 1000 different noise 
realizations. 



ionization gets enhanced by more than one order of mag- 
nitude, as compared to the case of the noise alone which 
is shown as the background featureless curve. Hence, the 
observed net enhancement can be attributed to a nonlin- 
ear quantum interaction between the coherent pulse and 
noise. 

As the noise level is further increased, we observe an 
enhancement of the periodic ionization profile by around 
three orders of magnitude as shown in Fig. 2(b). How- 
ever, the increase in noise level also causes the back- 
ground structureless stochastic ionization curve to rise 
monotonically. For strong noise (Fig. 2(c)), these peri- 
odic structures tend to wash out and the process is ef- 
fectively controlled by the noise. Hence one expects, the 
existence of an intermediate noise level where the nonlin- 
ear ionization is optimally enhanced. 

The net enhancement of the atomic ionization due to 
interplay between the laser pulse and the noise can be 
characterized by the enhancement factor 21], 



Pl+n — Pq 



(10) 



We have seen that the atomic photoionization due to 
an intense femtosecond laser pulse is a highly nonlinear 
quantum phenomenon, and in particular, the ionization 
response collapses for a "weak" laser pulse (see inset of 
Fig. 1). Motivated by the quantum SR effect, we wish 
to explore if the noise can recover the strong periodic 
ionization flux for the weak laser pulse. To answer this 
question, in Fig. 2(a) we show the average ionization 
flux when a small noise of amplitude, y/~D = 0.00024, is 
added to the previously weak laser pulse (F — 0.02). 
Note that the atomic excitation time by the laser and 
the noise here are identical. One can see that for such a 
feeble noise amplitude, the periodic structure in atomic 



with Po = Pi + P n . Although this is different compared 
to the quantifiers commonly used [l|,l2|], r\ is more suitable 
for our case. One can verify that a zero value of r\ corre- 
sponds to the case when either the laser pulse (Pi 3> P n ) 
or the noise (Pi <C P n ) dominates. In Fig. 4, we have 
plotted the enhancement factor r) versus the noise ampli- 
tude \[T). It exhibits a sharp rise, followed by a maxi- 
mum at a certain value of the noise (point B), and then a 
gradual fall off. It is worth mentioning that only a mod- 
est noise-to-laser ratio ( ■J D op t / Fq — 0.075) is required 
to reach the optimum enhancement (here rj max = 36). 

Before investigating other properties of the enhance- 
ment effect, it is worth making three remarks. First, 
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although the enhancement curve bears striking resem- 
blance to the typical SR curve, it is not SR effect where 
the matching of the time scales between coherent and in- 
coherent driving exists. The underlying gain mechanism 
here is completely different, as we shall see later. One can 
perhaps call this as a generalized quantum SR for such 
atomic systems, in the sense of the existence of an op- 
timum noise level. Second, the location of the optimum 
enhancement is governed by an empirical condition, when 
the strengths of the laser pulse and noise are comparable, 
in terms of the ionization flux produced by their individ- 
ual action P; ~ P n . This can be verified in Fig. 3 for the 
enhancement curve shown in Fig. 4. Third, due to the 
presence of the random noise term in the Hamiltonian, 
the optimal solution is only statistically unique. We have 
computed the standard deviation of the enhancement fac- 
tor rj using 1000 realizations, a = y~< rj 2 > — < r\ > 2 . 
The corresponding a values are shown by error bars on 
the rj curve in Fig. 4. 



^. Enhancement curves for a variety of laser pulses 

In this subsection, we study the sensitivity of the 
stochastic enhancement curves on the laser pulse. In 
particular, we study the role of two parameters: (i) the 
pulse duration T pi and (ii) the peak pulse amplitude Fq. 
In Fig. 5, we have plotted enhancement curves versus the 
noise amplitude for pulses of fixed amplitude (F = 0.02) 
but of varying durations from 5 to 30 optical cycles. One 
can clearly see that, the enhancement features (particu- 
larly the location and strength of the optima) are robust 
for pulses ranging from ultrashort few cycle duration to 
quite long ones. 

Although we do not show, we have also verified that 
variation in the carrier envelop phase 5 of the laser pulse 
F(t) [see Eq. (3)] does not modify the enhancement ef- 
fect. This can be expected due to the presence of the 
noise term, by which any effect of 6 is averaged out. 
Furthermore, we have also observed similar enhance- 
ment curves for other forms of the pulse envelop, such 
as /(*) = sin 2 (7rf/Tp). 

To investigate the dependence of rj on the laser pulse 
amplitude Fq, we choose some moderate noise amplitude 
value, for example, \f~D — 0.0015. For this fixed y 7 !?, we 
now increase the peak pulse amplitude Fq of the 20 cy- 
cles pulse from zero to a large value, and calculate the IP 
for each value of Fq. The obtained probabilities for dif- 
ferent cases are plotted in Fig. 6(a), which are then used 
to compute the enhancement factor rj in Fig. 6(b). Here, 
r\ also exhibits a nonmonotonic feature versus the laser 
peak amplitude, thus suggesting a range of Fo where the 
addition of noise can be useful. This dependency is intu- 
itively explained. Since for weak laser pulse the process is 
dominated by the noise and the rj collapses. On the other 
hand, if the laser peak amplitude is too strong (compa- 
rable to over-the-barrier ionization threshold), the pulse 
can ionize the atom by itself, and the noise has no role 
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FIG. 5: The enhancement factor versus the noise amplitude 
for four different values of pulse durations T p : 5, 10, 20, 30 
optical period To. Here Fo = 0.02, and other parameters are 
the same as in previous figures. 



to play. Thus, it is for the intermediate values of the 
noise and laser pulse amplitudes, where this nonlinear 
enhancement mechanism can be significant. From Figs. 4 
and 6, one can conclude that in order to maximize the 
net ionization yield, a particular pair of Fq and \f~D is 
required. 



B. Employing chaotic light instead of the white 
noise 

1. Generation and characterization of chaotic light 

To experimentally observe this effect, the most chal- 
lenging task is the generation of intense white noise. For 
instance, if one considers the thermal radiation from a 
blackbody (such as the sun) as a possible source of the 
white noise, its noise intensity falls short by many or- 
ders of magnitude [2(| , compared to the one required for 
the optimum of Fig. 4. We thus look into alternatives 
for generating a noise-like waveform. One possibility is 
to employ modern pulse shaping techniques, whereby one 
can design waveforms of almost arbitrary shapes [23, HI] . 
To realize such a chaotic light, we choose a large number 
of frequency modes, N , in a finite but broad bandwidth 
Auj. These modes can be, for example, different Fourier 
components of an ultrashort laser pulse. The total elec- 
tric field Z(t) is a sum of these N individual modes as 



ry n 

n=l 



, sm(oj n t + <f> n 



(11) 



4> n denote the angular frequency, phase of 
mode, respectively; and F rms is the root-mean- 
square amplitude of Z(t). Note that here we consider 



where 

„th 
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FIG. 6: The dependence of enhancement on Fq for a fixed 
noise amplitude \/T> = 0.0015. (a) Ionization probabilities 
versus Fo for laser only Pi, and for laser with noise P;+ n . 
The dashed line shows the limiting value of P n . (b) The 
enhancement vs Fo curve obtained from the probability curves 
shown in (a). 
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FIG. 7: Enhancement induced by a broadband chaotic light. 
The peak amplitude and frequency of the 20 cycle laser pulse 
are Fq = 0.02 and u = 0.057, respectively. The bandwidth 
of chaotic light is Auj = 0.75 with central frequency luq = 
0.375. Inset: power spectral density (PSD) of the chaotic 
light compared to the one for the white noise. 



light such as its dependence on the pulse duration T p 
and pulse amplitude Fo. 

This observation not only suggests the possibility of 
observing the effect using a finite but broadband chaotic 
light, but also raises the question concerning the relevant 
frequency components in the chaotic light spectrum, to 
be discussed next. 



Spectral and temporal analysis of the gain 
mechanism 



these frequency modes to oscillate independently with 
their phases </>„ assuming random values relative to each 
other. In this particular case of phase-randomized co- 
herent modes, the total field Z(t) at any point will be 
noise-like, fluctuating in intensity due to the interference 
between modes. Inset in Fig. 7 shows an example of such 
a chaotic light spectrum for N — 1024 in a chosen band- 
width (BW) of 0.75 (corresponding to a 32 attosecond 
pulse) [3fJ. Such a construct tends to the white noise, 
in the limit of Alu, N — > oo. In the following subsection, 
we consider a simultaneous application of the weak laser 
pulse and chaotic light (instead of the white noise) and 
see if the enhancement effect can be preserved. 



2. Photoionization by the chaotic light 

When we replace white noise by the chaotic light 
(BW=0.75 and N = 1024), one can see in Fig. 7 that 
most of the features of the enhancement phenomenon re- 
main intact. In particular, the intensity and location of 
the optimum is very close to the one obtained for the 
white noise case in Fig. 4. We also recover other fea- 
tures of the enhancement mechanism with the chaotic 



1. Frequency-resolved gain profile 

In this subsection, we will analyze the mechanism of 
stochastic enhancement in both, the frequency domain 
and time domain. In particular, we aim to identify fre- 
quency components in the broad spectrum of noise (or 
chaotic light) which are the crucial ones to provide the 
gain. For this purpose, we compute a frequency-resolved 
atomic gain (FRAG) profile using a pump-probe type of 
setting as described below. 

It is well known that when an atom interacts linearly 
with a weak external field, the energy absorption takes 
place at its resonant frequencies. However, due to its 
interaction with a strong laser pulse, the unperturbed 
atomic states are significantly modified leading to a com- 
pletely different frequency response of the driven atom. 
Indeed, such a modified spectral response is the relevant 
quantity when the atom is also subjected to noise. So, 
how can one measure precisely the frequency resolved 
gain G{v) offered by the atom? One possible way to 
observe FRAG is to consider the previously employed 
laser pulse (u = 0.057, Fo = 0.02) as a pump pulse and 
replace noise by a tunable monochromatic probe pulse, 
F p (t) = f(t)F p sin(cjpt). The probe amplitude F p is much 
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FIG. 8: Frequency-resolved gain G(uj p ) of the atom which is 
driven by a 10 cycle laser pulse of Fo = 0.02 at uj = 0.057. The 
atomic gain is probed by a simultaneous application of a weak 
probe beam F p = 0.0002 of tunable frequency ui p . Both, the 
depletion of the ground state population (blue) and the net 
absorption of energy (red) due to the probe are plotted. The 
gain of driven atom is also compared with the bare atomic 
gain (black). 
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FIG. 9: Enhancement curves due to chaotic light spectrum 
perforated by holes at the first three resonance frequencies. 
Difference curves correspond to the increasing hole widths of 
w = 0.0, 0.013, 0.03, 0.1 and 0.15. 



weaker than the driving laser pulse [F p /Fq = 0.01], such 
that it does not significantly alter the quasi energy lev- 
els, but can drive transitions between them. For an atom 
prepared in its ground state, a FRAG profile is obtained 
by measuring cither the depletion of the ground state 
population or the net energy absorption, as a function of 
the probe frequency uj p . 

Such a gain profile G(oj p ) is shown in Fig. 8 for the 
atom driven by a 10 cycle long laser pulse of amplitude 
F = 0.02. Although the FRAG shows lots of structure, 
the dominant peaks appear around the first atomic tran- 
sition frequency. A closer inspection of Fig. 8 reveals that 
there is a dip at the unperturbed resonance frequency 
ujqi = 0.267. The main peak is indeed shifted beyond the 
linear Stark shift for our model. These shifts and broad- 
ening of the atomic resonances are caused by the strong 
laser pulse, since the field can drive the electron signifi- 
cantly away from the nucleus leading to strong nonlinear 
perturbation to its bound states. Such a FRAG curve 
provides a fingerprint of the atomic gain under the strong 
laser pulse. The frequency bands that correspond to the 
peaks in the obtained gain curve are indeed the most 
useful ones to obtain the enhancement. In the following, 
we will test the validity of this statement by designing a 
chaotic light where significant resonance frequencies are 
filtered from its spectrum. 



hancement is lost, if the hole width w is below 0.013 a.u., 
which already includes the linear ac stark shift of the 
atomic states. However, by increasing the hole width 
such that no frequency component exists in the noise 
spectrum where the FRAG has dominant peaks leads 
to a collapse of the enhancement mechanism, as shown 
in Fig. 9. This observation validates the importance of 
FRAG structure in identifying the useful spectral bands 
in noise. It suggests that for the simultaneous presence 
of the laser pulse and noise, non-resonant frequency com- 
ponents are the dominant ones. 

It is worth making few remarks here. First, the FRAG 
is a property of the atom interacting with a particular 
strong laser pulse. Thus, the detailed features of the 
gain curve depends on both the atomic system and on the 
laser pulse parameters. But this doesn't affect the general 
conclusion drawn from such a curve. Second, it is also 
possible to obtain the enhancement using a monochro- 
matic beam instead of the chaotic broadband light, if 
its frequency is properly tuned to the new "resonances" . 
But, the advantage of using a broadband source is that 
the enhancement becomes independent of both, the par- 
ticular atom and the FRAG structure for different pulse 
parameters. 

3. Role of relative delay between signal and noise 



2. Chaotic light with missing resonant frequencies 

To show that the useful frequencies are not simply the 
resonant frequencies of the atom, we have designed a 
chaotic light spectrum perforated by digging holes in its 
spectral density around the first few atomic transition 
frequencies. As can be seen from Fig. 9, almost no en- 



Up to now we have considered the case of a perfect syn- 
chronization, i.e., a simultaneous application of the laser 
pulse and noise to observe the enhancement. We now 
wish to relax this synchronization constrain between the 
laser pulse and noise to test if the enhancement still ex- 
ists. Such a scenario would not only help the experimen- 
tal search of similar effects, but also provides an alterna- 
tive aspect of the enhancement mechanism as compared 
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FIG. 10: Enhancement factor r\ versus the noise amplitude 
\[D for various cases, application of noise first and then the 
laser pulse (empty squares) , application of laser pulse first and 
then noise (filled square), simultaneous application of laser 
and noise (circles). Insets depicts schematically temporal sig- 
nal applied to the atom for corresponding cases. 

to the one mentioned above. 

There are two possible ways to expose the atom to a 
laser pulse and noise sequentially: (i) the atom first in- 
teracts with noise only and then we apply the laser pulse, 
or conversely, (ii) the laser pulse is applied first and then 
the noise is applied. Note that for both cases, there is no 
direct interplay between laser and noise. The enhance- 
ment factor 77, which is defined as before, is shown in 
Fig. 10 for both the cases. One can clearly see that for 
the case-(i) r\ curve looks very similar to the one for simul- 
taneous action of laser and noise. But for the case-(ii), 
the enhancement curve collapses. 

Although we obtain almost identical enhancement 
curves for cases of, first-noise then laser pulse and si- 
multaneous action of noise with laser pulse, the gain pro- 
viding frequency components are fundamentally different 
in each case. The gain curves for both cases are given in 
Fig. 8. For the sequential application of noise and laser, 
the atomic gain is basically at the resonant frequencies 
of the unperturbed atom. Thus the resonant frequencies 
(particularly the first few) are the most fertile ones in the 
noise spectrum. In this case, a two-step picture of the 
enhancement mechanism applies, where the atom first 
absorbs energy from the noise leading to an exponential 
population distribution, and the laser causes ionization 
from "noise-heated" atom in a second step. 



IV. SUMMARY AND CONCLUSION 

We have investigated photoionization of a hydrogen 
atom which is subjected to both, a MIR femtosecond 
laser pulse and white Gaussian noise. Due to the in- 
herent nonlinearity of the ionization process, a form of 
quantum stochastic resonance-like behavior has been ob- 



served. This quantum SR leads to a dramatic enhance- 
ment (by several orders of magnitude) in the nonlinear 
ionization when a specific but small amount of white 
noise is added to the weak few cycle laser pulse. We have 
further shown the signatures of the enhancement effect 
for different types of the laser pulses from a few cycles to 
few tens of cycles duration. Moreover, if the noise am- 
plitude is kept fixed to some level, and the peak pulse 
amplitude is varied, again a curve with a specific maxi- 
mum is obtained for the enhancement parameter. These 
results suggest the existence of an optimum combination 
of the laser pulse and noise, if one is interested in optimiz- 
ing the relative ionization enhancement. The same effect 
is also achieved if one uses realizable broadband chaotic 
light instead of white noise. We emphasize that the ef- 
fect is robust with respect to a range of experimentally 
accessible parameters such as the pulse duration. 

The enhancement mechanism is analyzed in the fre- 
quency domain by measuring the frequency-resolved gain 
profile of the atom under a strong laser pulse, employing 
a pump-probe type of setting. The gain providing fre- 
quencies are significantly modified from the unperturbed 
atomic resonances, suggesting the non-resonant nature 
of the noise absorption. However, if we introduce a rel- 
ative time-delay between the laser pulse and noise the 
enhancement is still present, provided the noise acts first 
on the atom. In this case, the useful frequencies in the 
noise spectrum are at the atomic resonances. 

Finally, analogous effects are also expected in other 
systems provided the following three conditions are full- 
filled: (i) The system has a single-well finite binding po- 
tential with multiple energy levels, (ii) it can be subjected 
to a (nonresonant) coherent optical driving, and (iii) it 
can be subjected to an incoherent perturbation. Since 
these conditions are sufficiently general, other systems 
(SQUIDs, molecules etc) might also display similar fea- 
tures [mil]. 
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APPENDIX 

In this appendix we briefly present an algorithm for 
the numerical simulation of the time-dependent stochas- 
tic Schrodinger equation, Eq. (6) in the text, where the 
Hamiltonian is given by Eq. (1). The properties of the 
white Gaussian noise £(t) are defined in Eq. (4)-(5). Our 
basic approach is to use the split-operator fast-Fourier 
transform (SOFFT) method due to Feit and Fleck (which 
is well known for the deterministic case) pol ] , and adapt it 
to the case when the Hamiltonian contains an additional 
stochastic term £(t). 
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Recalling briefly that if there were no random term, 
the solution of Eq. (6) can be obtained by defining the 
standard propagator U(to, t), which when applied on ini- 
tial state wavefunction \^o) propagates it forward in 
time. The usual approach to compute the solution of 
the Schrodinger equation is to discretize the total propa- 
gation time into N small steps of equal intervals At. The 
resulting exact short-time propagator can be written as, 



Uz(At) = exp 



t+At 



(H det (x,t) + x£(t))dt 



(A-l) 

Here, Hd e t(x, t) — p 2 /2 + V(x) + xF(t), is the determinis- 
tic atomic Hamiltonian including the laser-atom interac- 
tion term. In order to incorporate the white noise term 
within the framework of SOFFT method, one can rewrite 
the propagator as, 



U 6 (At) = U (At) exp 



t+At 



x£{t)dt 



(A-2) 



where, U(At) = exp (— iHdet(%, t)At), denotes the deter- 
ministic part of the propagator. The stochastic integral 
in the exponential can be interpreted in the Stratonovitch 
sense |23j using the properties of the white Gaussian noise 
as follows, 



/t+At 
i(t)dt = V2DAtX t 



(A-3) 



where X t is a random number having Gaussian distribu- 
tion and of unit variance. This makes the propagator a 
stochastic operator. Note that even in the presence of 



the noise term the operator is unitary, i.e., it preserves 
the norm of the wavefunction. 

One can approximate the exact propagator given by 
Eq. (A-2) following a three-step splitting leading to the 
following expression, 



U$(At) = exp ( -*yAi 



exp (-iV((x, t)At)) exp ( -«y Ai 



(A-4) 



with, V s (x, t) = V(x) + xF(t) + y/2D/AtX t . Note that 
the effect of noise is simulated by inserting at every time 
step a random number X t whose statistical properties 
are described above. The right hand side of Eq. (A-4) is 
thus equivalent to the free propagation over a half time 
increment At/2, a random phase change from the action 
of potential V^(x,t) over the whole time At, and an ad- 
ditional free particle propagation over At/2. 

This operator splitting is correct up to second order in 
the time step At for the noise- free part, but due to the 
stochastic integration it is accurate up to only first-order 
for the stochastic part. In the actual calculation At is 
chosen sufficiently small, such that a further reduction 
in its value does not alter the accuracy of the physical 
results. For a given realization of the random number se- 
quence entering in the propagator via X t , one generates 
a quantum "trajectory" for the wavefunction. To extract 
the physical observable, an ensemble average of the de- 
sired quantity over a large number of noise realizations 
is needed. Other simulation parameters such as the grid 
size and the grid resolutions should be taken as described 
in the literature [241 ]. 
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